 <html>
 <head>
 <title>k_tpsalie_analysis.f90</title>                                          
                                                                                
                                           
 </head>
 <body>
MODULE TPSALIE_ANALYSIS <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  USE TPSALIE <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
  IMPLICIT NONE <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
  PUBLIC <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
  PRIVATE ALLOCVECRES,ALLOCPBRES,ALLOCONELIE,ALLOCDF,ALLOCFD,ALLOCNORMAL <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
  PRIVATE NORMALMAP,MAPREVDF,REVDFMAP,MAPDF,DFMAP,RESPB,PBRES,RESOVEC,VECRESO <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
  PRIVATE ONEEXPMAP,KILLNORMAL,KILLDF,KILLREVDF,KILLPBRES,KILLVECRES,KILLONELIE <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
  PRIVATE DAPRINTVECRES,DAPRINTPBRES,DAPRINTDF,DAPRINTREVDF <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
  PRIVATE DF_MAP,FD_MAP,MAP_DF,MAP_FD,ONE_MAP,MAP_ONE,MAPNORMAL <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
  PRIVATE  RESTA,TARES,ALLOCTARES,KILLTARES <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
  PRIVATE ALLOCGEN,KILLGEN <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
  PRIVATE EQUALGENMAP,EQUALMAPGEN,PUSHGEN <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
  INTEGER,PRIVATE::NO,ND,ND2,NP,NDPT,NV <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  LOGICAL(LP),PRIVATE::OLD <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
  LOGICAL(LP) IMAXFLAG <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         
  PRIVATE NRESO <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="INIT_TPSALIE"></a>
<font color="#0000FF">  INTERFACE INIT_TPSALIE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#INIT_MAP">     MODULE PROCEDURE INIT_MAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#INIT_TPSA">     MODULE PROCEDURE INIT_TPSA</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="="></a>
<font color="#0000FF">  INTERFACE ASSIGNMENT (=) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#NORMALMAP">     MODULE PROCEDURE NORMALMAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAPNORMAL">     MODULE PROCEDURE MAPNORMAL</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAPREVDF">     MODULE PROCEDURE MAPREVDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#REVDFMAP">     MODULE PROCEDURE REVDFMAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAPDF">     MODULE PROCEDURE MAPDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DFMAP">     MODULE PROCEDURE DFMAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#RESPB">     MODULE PROCEDURE RESPB</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#PBRES">     MODULE PROCEDURE PBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#RESTA">     MODULE PROCEDURE RESTA</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#TARES">     MODULE PROCEDURE TARES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#RESOVEC">     MODULE PROCEDURE RESOVEC</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#VECRESO">     MODULE PROCEDURE VECRESO</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAPONEEXP">     MODULE PROCEDURE MAPONEEXP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ONEEXPMAP">     MODULE PROCEDURE ONEEXPMAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#EQUALMAPGEN">     MODULE PROCEDURE EQUALMAPGEN</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#EQUALGENMAP">     MODULE PROCEDURE EQUALGENMAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font color="#0000FF">  END  INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="*"></a>
<font color="#0000FF">  INTERFACE OPERATOR (*) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DF_MAP">     MODULE PROCEDURE DF_MAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#FD_MAP">     MODULE PROCEDURE FD_MAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAP_DF">     MODULE PROCEDURE MAP_DF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAP_FD">     MODULE PROCEDURE MAP_FD</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ONE_MAP">     MODULE PROCEDURE ONE_MAP</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#MAP_ONE">     MODULE PROCEDURE MAP_ONE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#PUSHGEN">     MODULE PROCEDURE PUSHGEN</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font color="#0000FF">  END  INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOC"></a>
<font color="#0000FF">  INTERFACE ALLOC </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCVECRES">     MODULE PROCEDURE ALLOCVECRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCPBRES">     MODULE PROCEDURE ALLOCPBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCONELIE">     MODULE PROCEDURE ALLOCONELIE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCDF">     MODULE PROCEDURE ALLOCDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCFD">     MODULE PROCEDURE ALLOCFD</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCNORMAL">     MODULE PROCEDURE ALLOCNORMAL</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCGEN">     MODULE PROCEDURE ALLOCGEN</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#ALLOCTARES">     MODULE PROCEDURE ALLOCTARES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                   
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILL"></a>
<font color="#0000FF">  INTERFACE KILL </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLONELIE">     MODULE PROCEDURE KILLONELIE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLNORMAL">     MODULE PROCEDURE KILLNORMAL</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLDF">     MODULE PROCEDURE KILLDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLREVDF">     MODULE PROCEDURE KILLREVDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLPBRES">     MODULE PROCEDURE KILLPBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLVECRES">     MODULE PROCEDURE KILLVECRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLGEN">     MODULE PROCEDURE KILLGEN</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#KILLTARES">     MODULE PROCEDURE KILLTARES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAPRINT"></a>
<font color="#0000FF">  INTERFACE DAPRINT </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTVECRES">     MODULE PROCEDURE DAPRINTVECRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTPBRES">     MODULE PROCEDURE DAPRINTPBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTDF">     MODULE PROCEDURE DAPRINTDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTREVDF">     MODULE PROCEDURE DAPRINTREVDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTONELIE">     MODULE PROCEDURE DAPRINTONELIE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                             
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="PRINT"></a>
<font color="#0000FF">  INTERFACE PRINT </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTVECRES">     MODULE PROCEDURE DAPRINTVECRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTPBRES">     MODULE PROCEDURE DAPRINTPBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTDF">     MODULE PROCEDURE DAPRINTDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTREVDF">     MODULE PROCEDURE DAPRINTREVDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAPRINTONELIE">     MODULE PROCEDURE DAPRINTONELIE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                             
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAINPUT"></a>
<font color="#0000FF">  INTERFACE DAINPUT </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADVECRES">     MODULE PROCEDURE DAREADVECRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADPBRES">     MODULE PROCEDURE DAREADPBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADDF">     MODULE PROCEDURE DAREADDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADREVDF">     MODULE PROCEDURE DAREADREVDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADONELIE">     MODULE PROCEDURE DAREADONELIE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="READ"></a>
<font color="#0000FF">  INTERFACE READ </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADVECRES">     MODULE PROCEDURE DAREADVECRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADPBRES">     MODULE PROCEDURE DAREADPBRES</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADDF">     MODULE PROCEDURE DAREADDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADREVDF">     MODULE PROCEDURE DAREADREVDF</a></font> <br>                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="#DAREADONELIE">     MODULE PROCEDURE DAREADONELIE</a></font> <br>                                                                                                                                                                                                                                                                                                                                                               
<font color="#0000FF">  END INTERFACE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
CONTAINS <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ONE_MAP"></a>
<font color="#37A50B">  FUNCTION ONE_MAP(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) ONE_MAP ,JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP), INTENT (IN) :: S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOCALMASTER=MASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECKDAMAP(S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSDAMAP(ONE_MAP)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ONE_MAP=JUNK*S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION ONE_MAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAP_ONE"></a>
<font color="#37A50B">  FUNCTION MAP_ONE(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) MAP_ONE,JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP), INTENT (IN) :: S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOCALMASTER=MASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECKDAMAP(S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSDAMAP(MAP_ONE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MAP_ONE=S2*JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION MAP_ONE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DF_MAP"></a>
<font color="#37A50B">  FUNCTION DF_MAP(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) DF_MAP ,JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP), INTENT (IN) :: S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOCALMASTER=MASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECKDAMAP(S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSDAMAP(DF_MAP)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DF_MAP=JUNK*S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION DF_MAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                            
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAP_DF"></a>
<font color="#37A50B">  FUNCTION MAP_DF(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) MAP_DF,JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP), INTENT (IN) :: S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOCALMASTER=MASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECKDAMAP(S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSDAMAP(MAP_DF)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MAP_DF=S2*JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION MAP_DF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                            
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="FD_MAP"></a>
<font color="#37A50B">  FUNCTION FD_MAP(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) FD_MAP ,JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP), INTENT (IN) :: S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOCALMASTER=MASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECKDAMAP(S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSDAMAP(FD_MAP)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    FD_MAP=JUNK*S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION FD_MAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                            
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAP_FD"></a>
<font color="#37A50B">  FUNCTION MAP_FD(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) MAP_FD ,JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP), INTENT (IN) :: S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOCALMASTER=MASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECKDAMAP(S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSDAMAP(MAP_FD)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MAP_FD=S2*JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=LOCALMASTER</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION MAP_FD </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                            
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="GO_TO_FIX_POINT"></a>
<font color="#FF0000">  SUBROUTINE  GO_TO_FIX_POINT(M,A1,A1I,NORD) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(INOUT)::M,A1,A1I</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER NORD</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL GOFIX(M%V%I,A1%V%I,A1I%V%I,NORD)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE GO_TO_FIX_POINT </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="NORMALMAP"></a>
<font color="#FF0000">  SUBROUTINE  NORMALMAP(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (NORMALFORM),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) ZERO_(NDIM2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOGICAL(LP) GLO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    GLO=GLOBAL_VERBOSE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ZERO_(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%NORMAL%LINEAR%V(1)%I==0)  CALL CRAP1("NORMALMAP 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LIELIB_PRINT(7)=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.S2%AUTO) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LIELIB_PRINT(7)=-1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL SETIDPR(S2%PLANE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL INPUTRES(S2%M,S2%NRES)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A%CONSTANT=JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=ZERO_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%NORD.LE.0.OR.S2%NORD.GT.NO) S2%NORD=NO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%JTUNE.LT.0) S2%JTUNE=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL MAPNORMF(JUNK%V%I,S2%A%NONLINEAR%V%I,S2%A%LINEAR%V%I,&</font> <br>                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;         & S2%A1%V%I,S2%NORMAL%LINEAR%V%I,S2%NORMAL%NONLINEAR%V%I,S2%NORD,S2%JTUNE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NORMAL%LINEAR=(S2%NORMAL%LINEAR).SUB.1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NORMAL%PB=S2%NORMAL%NONLINEAR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A%PB=S2%A%NONLINEAR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DHDJFLO(S2%NORMAL%NONLINEAR%V%I,S2%DHDJ%V%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL GETTURA(S2%TUNE,S2%DAMPING)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A_T=S2%A</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A_T=S2%A1*S2%A_T</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A_T=ZERO_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    GLOBAL_VERBOSE=GLO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE NORMALMAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAPNORMAL"></a>
<font color="#FF0000">  SUBROUTINE  MAPNORMAL(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (NORMALFORM),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) JUNK,ID</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK,ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ID=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ID=TEXP(S2%NORMAL%NONLINEAR,ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S2%A1*S2%A</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=JUNK**(-1)*ID*JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK,ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE MAPNORMAL </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ONEEXPMAP"></a>
<font color="#FF0000">  SUBROUTINE  ONEEXPMAP(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%VECTOR%V(1)%I==0) CALL CRAP1("ONEEXPMAP 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%EPS==0.0_DP) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       S2%EPS=1E3_DP*FULL_ABS(S1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL FLOFACG(JUNK%V%I,S2%VECTOR%V%I,S2%EPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       S2%EPS=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ELSE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL FLOFACG(JUNK%V%I,S2%VECTOR%V%I,S2%EPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB=S2%VECTOR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font color="#FF0000">  END SUBROUTINE ONEEXPMAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAPONEEXP"></a>
<font color="#FF0000">  SUBROUTINE  MAPONEEXP(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=TEXP(S2%VECTOR,JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font color="#FF0000">  END SUBROUTINE MAPONEEXP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="REVDFMAP"></a>
<font color="#FF0000">  SUBROUTINE  REVDFMAP(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) ZERO_(NDIM2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ZERO_(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%LINEAR%V(1)%I==0)  CALL CRAP1("REVDFMAP 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%CONSTANT=JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=ZERO_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=JUNK**(-1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL FLOFAC(JUNK%V%I,S2%LINEAR%V%I,S2%NONLINEAR%V%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%LINEAR=S2%LINEAR**(-1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DACMUD(S2%NONLINEAR%V%I,-1.0_DP,S2%NONLINEAR%V%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB=S2%NONLINEAR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE REVDFMAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DFMAP"></a>
<font color="#FF0000">  SUBROUTINE  DFMAP(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) ZERO_(NDIM2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ZERO_(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%LINEAR%V(1)%I==0)  CALL CRAP1("DFMAP 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%CONSTANT=JUNK</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    JUNK=ZERO_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL FLOFAC(JUNK%V%I,S2%LINEAR%V%I,S2%NONLINEAR%V%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB=S2%NONLINEAR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font color="#FF0000">  END SUBROUTINE DFMAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="RESOVEC"></a>
<font color="#FF0000">  SUBROUTINE  RESOVEC(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECFIELD),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECRESONANCE),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S1%COS%V(1)%I==0) CALL CRAP1("RESOVEC 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CTORFLO(S2%V%I,S1%COS%V%I,S1%SIN%V%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%IFAC=S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%COS%IFAC=S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%SIN%IFAC=S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE RESOVEC </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="VECRESO"></a>
<font color="#FF0000">  SUBROUTINE  VECRESO(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECFIELD),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECRESONANCE),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%V(1)%I==0) CALL CRAP1("VECRESO 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL RTOCFLO(S1%COS%V%I,S1%SIN%V%I,S2%V%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%IFAC=S1%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font color="#FF0000">  END SUBROUTINE VECRESO </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="RESPB"></a>
<font color="#FF0000">  SUBROUTINE  RESPB(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBFIELD),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBRESONANCE),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S1%COS%H%I==0) CALL CRAP1("RESPB 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CTOR(S2%H%I,S1%COS%H%I,S1%SIN%H%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%IFAC=S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%COS%IFAC=S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%SIN%IFAC=S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE RESPB </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="PBRES"></a>
<font color="#FF0000">  SUBROUTINE  PBRES(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBFIELD),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBRESONANCE),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%H%I==0) CALL CRAP1("RESPB 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL RTOC(S1%COS%H%I,S1%SIN%H%I,S2%H%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%IFAC=S1%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font color="#FF0000">  END SUBROUTINE PBRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="RESTA"></a>
<font color="#FF0000">  SUBROUTINE  RESTA(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (TAYLOR),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (TAYLORRESONANCE),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S1%COS%I==0)  CALL CRAP1("RESTA 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CTOR(S2%I,S1%COS%I,S1%SIN%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font color="#FF0000">  END SUBROUTINE RESTA </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="TARES"></a>
<font color="#FF0000">  SUBROUTINE  TARES(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (TAYLOR),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (TAYLORRESONANCE),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%I==0)  CALL CRAP1("TARES 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL RTOC(S1%COS%I,S1%SIN%I,S2%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font color="#FF0000">  END SUBROUTINE TARES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAPDF"></a>
<font color="#FF0000">  SUBROUTINE  MAPDF(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) ID</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER NO1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NO1=NO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S1%V(1)%I==0) CALL CRAP1("MAPDF 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ID=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ID=TEXPDF( S2%NONLINEAR, ID,2,NO1,1.0_DP,1 )</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=ID*S2%LINEAR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=S2%CONSTANT</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font color="#FF0000">  END SUBROUTINE MAPDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="MAPREVDF"></a>
<font color="#FF0000">  SUBROUTINE  MAPREVDF(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP) ID</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER NO1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NO1=NO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S1%V(1)%I==0) CALL CRAP1("MAPREVDF 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ID=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ID=TEXPDF( S2%NONLINEAR, ID,2,NO1,1.0_DP,-1 )</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=S2%LINEAR*ID</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=S2%CONSTANT</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(ID)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font color="#FF0000">  END SUBROUTINE MAPREVDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCGEN"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCGEN(S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (GENFIELD),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I,J</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(S1%H)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%EPS=1E-6_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%IFAC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%LINEAR_IN=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%IMAX=1000</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%CONSTANT(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMAXFLAG=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%NO_CUT=NO+1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO J=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          CALL ALLOC(S1%D(I,J))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL ALLOC(S1%LINEAR%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL ALLOC(S1%LINEART%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL ALLOC(S1%M%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL ALLOC(S1%MT%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE ALLOCGEN </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLGEN"></a>
<font color="#FF0000">  SUBROUTINE  KILLGEN(S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (GENFIELD),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I,J</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%IMAX=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1%NO_CUT=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO J=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          CALL KILL(S1%D(I,J))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(S1%LINEART%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(S1%LINEAR%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(S1%MT%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(S1%M%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DADAL1(S1%H%I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font color="#FF0000">  END SUBROUTINE KILLGEN </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="EQUALGENMAP"></a>
<font color="#FF0000">  SUBROUTINE  EQUALGENMAP(S2,S1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (GENFIELD),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP)  W</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE(TAYLOR) T</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE(ONELIEEXPONENT) ONE_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) ZERO_(NDIM2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I,J,JN(LNV),K</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(W)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(T)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(ONE_)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%CONSTANT=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,NDIM2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ZERO_(I)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,LNV</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       JN(I)=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       JN(2*I)=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W=S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W=ZERO_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%LINEAR=W.SUB.1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%LINEAR_IN) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       S2%LINEAR=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W=W*S2%LINEAR**(-1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL MAKETREE(S2%LINEAR%V(I),S2%LINEART%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%IFAC>1) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ONE_=W</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO K=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          ONE_%VECTOR%V(K)=ONE_%VECTOR%V(K)/S2%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       W=ONE_</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%H%I==0) CALL CRAP1("EQUALGENMAP 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S2%M%V(1)%I==0) CALL CRAP1("EQUALGENMAP 2")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ETPIN(W%V%I,S2%M%V%I,JN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL INTD(S2%M%V%I,S2%H%I,1.0_DP)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%M=S2%M.CUT.S2%NO_CUT</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO J=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          T=  (S2%M%V(2*I)).D.(2*J)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          CALL MAKETREE(T,S2%D(I,J))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL MAKETREE(S2%M%V(I),S2%MT%V(I))</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(ONE_)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(T)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(W)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE EQUALGENMAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="EQUALMAPGEN"></a>
<font color="#FF0000">  SUBROUTINE  EQUALMAPGEN(S1,S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (GENFIELD),INTENT(IN)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DAMAP)  W</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I,JN(LNV)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL CHECK_SNAKE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(W)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,LNV</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       JN(I)=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       JN(2*I)=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(S1%V(1)%I==0)CALL CRAP1("EQUALMAPGEN 1")</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DIFD(S2%H%I,W%V%I,1.0_DP)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ETPIN(W%V%I,S1%V%I,JN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=S1*S2%LINEAR</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S1=S2%CONSTANT</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILL(W)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE EQUALMAPGEN </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="PUSHGEN"></a>
<font color="#37A50B">  FUNCTION PUSHGEN( S1, S2 ) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (GENFIELD), INTENT (IN) :: S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP), INTENT(IN),DIMENSION(:)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) PUSHGEN(LNV),S2T(NV)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) JUNK(LNV),JUNK2(LNV),ET(NDIM),DE(NDIM),MT(NDIM,NDIM),MI(NDIM,NDIM)</font> <br>                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP) EB,E</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOGICAL(LP) MORE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I,J,K,IMAX,IER,IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.C_%STABLE_DA) RETURN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(.NOT.IMAXFLAG) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       JUNK(:)=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          JUNK(I)= PUSH1POL(S1%LINEART%V(I),S2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO I=ND2+1,NV</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          JUNK(I)=S2(I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       IMAX=S1%IMAX</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO IFAC=1,S1%IFAC</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          MORE=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          EB=1.1111E5_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          DO I=1,NV</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             S2T(I)= JUNK(I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          DO I=1,IMAX</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             DE(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             DO K=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                DE(K)=S2T(2*K)-PUSH1POL(S1%MT%V(2*K),JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             DO J=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                DO K=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   MT(J,K)=PUSH1POL(S1%D(J,K),JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             CALL MATINV(MT,MI,ND,NDIM,IER)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ET(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             DO J=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                DO K=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   ET(J)=MI(J,K)*DE(K)+ET(J)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             E=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             DO K=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                JUNK(2*K)=JUNK(2*K)+ET(K)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                E=ABS(ET(K))+E</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             IF(MORE) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                IF(E>S1%EPS) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   EB=E</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                ELSE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   EB=E</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   MORE=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ELSE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                IF(E>=EB) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   GOTO 100</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                ELSE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                   EB=E</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;                ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          IMAXFLAG=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          GOTO 1000</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;100       CONTINUE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          DO K=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             JUNK2(2*K-1)=PUSH1POL(S1%MT%V(2*K-1),JUNK)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          DO K=1,ND</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;             JUNK(2*K-1)=JUNK2(2*K-1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO K=1,NV</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          PUSHGEN(K)=JUNK(K)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;1000 CONTINUE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#37A50B">  END FUNCTION PUSHGEN </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCPBRES"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCPBRES(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBRESONANCE),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(S2%COS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(S2%SIN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font color="#FF0000">  END SUBROUTINE ALLOCPBRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCTARES"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCTARES(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (TAYLORRESONANCE),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(S2%COS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(S2%SIN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font color="#FF0000">  END SUBROUTINE ALLOCTARES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCVECRES"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCVECRES(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECRESONANCE),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%COS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%SIN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font color="#FF0000">  END SUBROUTINE ALLOCVECRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCDF"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCDF(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NONLINEAR%IFAC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB%IFAC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%CONSTANT(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font color="#FF0000">  END SUBROUTINE ALLOCDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCONELIE"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCONELIE(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT)  S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%VECTOR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%EPS=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font color="#FF0000">  END SUBROUTINE ALLOCONELIE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCNORMAL"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCNORMAL(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (NORMALFORM),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I,J</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%A_T)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%A%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%A%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%A%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%NORMAL%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%NORMAL%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%NORMAL%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NORMAL%CONSTANT(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A%CONSTANT(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NORMAL%NONLINEAR%IFAC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NORMAL%PB%IFAC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A%NONLINEAR%IFAC=-1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%A%PB%IFAC=-1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NORD=NO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%AUTO=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NRES=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%JTUNE=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,NDIM</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       S2%TUNE(I)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       S2%DAMPING(I)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       S2%PLANE(I)=2*I-1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       DO J=1,NRESO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;          S2%M(I,J)=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%A1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%DHDJ)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font color="#FF0000">  END SUBROUTINE ALLOCNORMAL </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="ALLOCFD"></a>
<font color="#FF0000">  SUBROUTINE  ALLOCFD(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOCTPSA(S2%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NONLINEAR%IFAC=-1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB%IFAC=-1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%CONSTANT(:)=0.0_DP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font color="#FF0000">  END SUBROUTINE ALLOCFD </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLNORMAL"></a>
<font color="#FF0000">  SUBROUTINE  KILLNORMAL(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (NORMALFORM),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%A_T)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%A%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%A%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%A%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%NORMAL%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%NORMAL%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%NORMAL%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%A1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%DHDJ)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font color="#FF0000">  END SUBROUTINE KILLNORMAL </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLONELIE"></a>
<font color="#FF0000">  SUBROUTINE  KILLONELIE(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT)  S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%VECTOR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font color="#FF0000">  END SUBROUTINE KILLONELIE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLPBRES"></a>
<font color="#FF0000">  SUBROUTINE  KILLPBRES(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBRESONANCE),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%COS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%SIN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font color="#FF0000">  END SUBROUTINE KILLPBRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLTARES"></a>
<font color="#FF0000">  SUBROUTINE  KILLTARES(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (TAYLORRESONANCE),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%COS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%SIN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font color="#FF0000">  END SUBROUTINE KILLTARES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLVECRES"></a>
<font color="#FF0000">  SUBROUTINE  KILLVECRES(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECRESONANCE),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%COS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%SIN)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font color="#FF0000">  END SUBROUTINE KILLVECRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLDF"></a>
<font color="#FF0000">  SUBROUTINE  KILLDF(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NONLINEAR%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE KILLDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                          
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILLREVDF"></a>
<font color="#FF0000">  SUBROUTINE  KILLREVDF(S2) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN),INTENT(INOUT)::S2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%LINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%NONLINEAR)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL KILLTPSA(S2%PB)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%NONLINEAR%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    S2%PB%IFAC=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE KILLREVDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAPRINTONELIE"></a>
<font color="#FF0000">  SUBROUTINE  DAPRINTONELIE(S1,MFILE,DEPS) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP),OPTIONAL,INTENT(INOUT)::DEPS</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    WRITE(MFILE,*) S1%EPS, " CONVERGENCE TEST NUMBER"</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%VECTOR,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%PB,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAPRINTONELIE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAREADONELIE"></a>
<font color="#FF0000">  SUBROUTINE  DAREADONELIE(S1,MFILE) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (ONELIEEXPONENT),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    READ(MFILE,*) S1%EPS</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%VECTOR,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%PB,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAREADONELIE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAPRINTVECRES"></a>
<font color="#FF0000">  SUBROUTINE  DAPRINTVECRES(S1,MFILE,DEPS) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECRESONANCE),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP),OPTIONAL,INTENT(INOUT)::DEPS</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%COS,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%SIN,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAPRINTVECRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                   
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAREADVECRES"></a>
<font color="#FF0000">  SUBROUTINE  DAREADVECRES(S1,MFILE) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (VECRESONANCE),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%COS,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%SIN,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAREADVECRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAPRINTDF"></a>
<font color="#FF0000">  SUBROUTINE  DAPRINTDF(S1,MFILE,DEPS) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP),OPTIONAL,INTENT(INOUT)::DEPS</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       WRITE(MFILE,*) S1%CONSTANT(I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%LINEAR,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%NONLINEAR,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAPRINTDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAREADDF"></a>
<font color="#FF0000">  SUBROUTINE  DAREADDF(S1,MFILE) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (DRAGTFINN),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       READ(MFILE,*) S1%CONSTANT(I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%LINEAR,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%NONLINEAR,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAREADDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAPRINTREVDF"></a>
<font color="#FF0000">  SUBROUTINE  DAPRINTREVDF(S1,MFILE,DEPS) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP),OPTIONAL,INTENT(INOUT)::DEPS</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       WRITE(MFILE,*) S1%CONSTANT(I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%LINEAR,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%NONLINEAR,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAPRINTREVDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAREADREVDF"></a>
<font color="#FF0000">  SUBROUTINE  DAREADREVDF(S1,MFILE) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (REVERSEDRAGTFINN),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER I</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    DO I=1,ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       READ(MFILE,*) S1%CONSTANT(I)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDDO</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%LINEAR,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%NONLINEAR,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAREADREVDF </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAPRINTPBRES"></a>
<font color="#FF0000">  SUBROUTINE  DAPRINTPBRES(S1,MFILE,DEPS) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBRESONANCE),INTENT(IN)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    REAL(DP),OPTIONAL,INTENT(INOUT)::DEPS</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%COS,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAPRINT(S1%SIN,MFILE,DEPS)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                        
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAPRINTPBRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DAREADPBRES"></a>
<font color="#FF0000">  SUBROUTINE  DAREADPBRES(S1,MFILE) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER,INTENT(IN)::MFILE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    TYPE (PBRESONANCE),INTENT(INOUT)::S1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%COS,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL DAINPUT(S1%SIN,MFILE)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DAREADPBRES </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="INIT_MAP"></a>
<font color="#FF0000">  SUBROUTINE INIT_MAP(NO1,ND1,NP1,NDPT1,LOG1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER NO1,ND1,NP1,NDPT1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOGICAL(LP) LOG1,PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NB_=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    PRESENT_TPSA=LINGYUN_YANG</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(LAST_TPSA==1) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DATERMINATE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ELSEIF(LAST_TPSA==2) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DATERMINATE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LINGYUN_YANG=PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL RESET_APERTURE_FLAG</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    OLD=LOG1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NO=NO1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ND=ND1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ND2=ND1*2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NP=NP1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NDPT=NDPT1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NV=ND2+NP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NEWPRINT=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL LIEINIT(NO1,NV,ND1,NDPT1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%NC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=(/" BERZ'S PACKAGE  "/)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%FC='(1((1X,A72),/))'</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL SET_IN_TPSA( NO,ND,ND2,NP,NDPT,NV,OLD)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL SET_IN_TPSALIE( NO,ND,ND2,NP,NDPT,NV,OLD)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%NC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=(/"          NO          ND         ND2          NP        NDPT          NV"/)</font> <br>                                                                                                                                                                                                                                                                                                                                                                     
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%FC='(1((1X,A72)))'</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%FI='(1X,6(6X,I6))'</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=(/NO,ND,ND2,NP,NDPT,NV/)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSIGN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NPARA_FPP=ND2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NPARA_ORIGINAL=NPARA_FPP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font color="#FF0000">  END SUBROUTINE INIT_MAP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                        
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="KILL_FPP"></a>
<font color="#FF0000">  SUBROUTINE KILL_FPP() </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOGICAL PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    PRESENT_TPSA=LINGYUN_YANG</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(LAST_TPSA==1) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DATERMINATE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DEALLOC_ALL</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ELSEIF(LAST_TPSA==2) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DATERMINATE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DEALLOC_ALL</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LINGYUN_YANG=PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font color="#FF0000">  END  SUBROUTINE KILL_FPP </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="INIT_TPSA"></a>
<font color="#FF0000">  SUBROUTINE INIT_TPSA(NO1,NP1,LOG1) </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    INTEGER NO1,ND1,NP1,NDPT1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOGICAL(LP) LOG1,PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NB_=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    PRESENT_TPSA=LINGYUN_YANG</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(LAST_TPSA==1) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DATERMINATE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ELSEIF(LAST_TPSA==2) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DATERMINATE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL KILL(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                    
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LINGYUN_YANG=PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    MASTER=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    OLD=LOG1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ND1=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NDPT1=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NO=NO1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ND=ND1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ND2=ND1*2</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NP=NP1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NDPT=NDPT1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                             
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NV=ND2+NP</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NEWPRINT=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL RESET_APERTURE_FLAG</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL LIEINIT(NO1,NV,ND1,NDPT1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%NC=1</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P=(/" BERZ'S PACKAGE  "/)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    W_P%FC='(1((1X,A72),/))'</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                               
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL SET_IN_TPSA( NO,ND,ND2,NP,NDPT,NV,OLD)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL SET_IN_TPSALIE( NO,ND,ND2,NP,NDPT,NV,OLD)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                         
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ASSIGN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(VARF1)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    CALL ALLOC(VARF2)</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                      
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NPARA_FPP=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                            
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    NPARA_ORIGINAL=0</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font color="#FF0000">  END SUBROUTINE INIT_TPSA </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 <a name="DATERMINATE"></a>
<font color="#FF0000">  SUBROUTINE DATERMINATE() </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IMPLICIT NONE</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                          
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LOGICAL PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                   
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    PRESENT_TPSA=LINGYUN_YANG</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    IF(LAST_TPSA==1) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.TRUE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                 
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DEASSIGN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ELSEIF(LAST_TPSA==2) THEN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       LINGYUN_YANG=.FALSE.</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       CALL DEASSIGN</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                       
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    ENDIF</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                  
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;    LINGYUN_YANG=PRESENT_TPSA</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                              
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font size= -1>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                           
<font color="#FF0000">  END SUBROUTINE DATERMINATE </font> <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                     
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
END MODULE TPSALIE_ANALYSIS <br>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    
<hr color="#CC0000" size="10" align="left" noshade>                                                                                                                                                                                                                                                                                                                                                                                                                                                                 
 </body>
 </html>
